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Abstract. The impact of various types of external potentials to the Peyrard- 



> 

■ Bishop DNA denaturation is investigated through statistical mechanics approach. 

The partition function is obtained using transfer integral method, and further the 
stretching of hydrogen bond is calculated using time independent perturbation method. 
It is shown that all types of external potentials accelerate the denaturation processes 
at lower temperature. In particular, it is argued that the Gaussian potential with 
infinitesimal width reproduces a constant force at one end of DNA sequence as already 
done in some previous works. 
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1. Introduction 

It has been known that the base-pairs of double-hehx break up and dissociate from each 
other to form two separated random coils when a solution of DNA macromolecules is 
heated up to 80°C. This phenomenon is referred as the DNA denaturation or the thermal 
DNA melting ([1] and references therein). On the other hand, at lower temperature 
before melting two strands of DNA can also be separated by applying an oppositely 
directed force on two strands at the DNA terminal point. This is known as the force- 
induced DNA melting [1]. Theoretically, the phenomena have been discussed in many 
works, for instance the work by Hankc et.al.on the denaturation of stretched DNA [2]. 
Also using the worm-like chain, the opened double-stranded DNA can be explained by a 
force exceeding certain critical value [3, 4]. On the other hand, the study of shear effect 
of pulling force shows that the shear unzipping of a heteropolymer would be similar to 
the unzipping in a tensile mode with sequence heterogeneity [5, 6]. 

It has been shown that a random external force would drastically change the phase 
diagram, while the ground state develops bubbles with various lengths as the random 
force fluctuation is increased. The fluctuating force denaturates the DNA by a gradual 
increase of bubble sizes. This suggests the possibility of opening up local bubbles 
in the selective regions without breaking the whole DNA [7]. The denaturation at 
lower temperature is realized by the enzymes, protein and so forth [8]. The effects 
are influenced by many factors such as the alkaline compound filling up a cell. Some 
works also investigated another effects, for instance the phase transition in a short DNA 
using the Peyrard- Bishop (PB) with additional delta function potential [9] . The paper 
argued that the denaturation is a localization-delocalization transition process. In the 
case without any external force, the denaturation is a second order transition, while 
under an external force it becomes the first-order one [1]. Following the same line, the 
present paper deals with external potentials interacting to DNA within the PB model, 
and investigates its effects to the denaturation process. 

This paper is organized as follows. First of all the PB model with various external 
potentials is brieffy introduced. It is then followed with the statistical mechanics 
formulation on the partition function of Hamiltonian under consideration. Before ending 
the paper with summary, the detail calculation of hydrogen bond stretching using first 
order time independent perturbation and its thermal behavior are presented in Sees. 3 
and 4. 

2. The Model 

Following the original PB model, the motion of DNA molecules is represented by 
transversal displacement, Un and Vn, corresponding to the base displacement from their 
equilibrium position along the direction of hydrogen bonds that is represented by the 
Morse potential to connect two bases in a pair [10]. 

The paper adopts the PB model with an additional external potential V which has 
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the same direction as Morse potential, 

n n 
+ f (^n - ^n-l)' + f (c^t^(— ) - l) V ViUr, - V^) . (1) 

Performing a transformation to the center of mass coordinate representing the in- 
phase and out-phase transversal motions, one can define — {un + Vn) / and 
Yn — {un — Vn)/V^- Then, it yields the PB Hamiltonian [10], 

n n 

+ f (Yn^i - Ynf + f (e-"/^^" - 1)' + V{Y^) , (2) 

where D and a are the depth and inverse width of the potential respectively. Pn = MXn, 
Pn = MYn and k is the spring constant. 

It should be remarked that some previous works suggested that the n— th nucleotide 
in one strand might preferably interacts with the {n ± h)—th nucleotides in another 
strand with h — A or 5 [11, 12, 13]. However for the sake of simphcity, throughout the 
paper let us consider the original PB model with taking into account only the nearest 
neighbor nucleotides of another strand. The extension to the cases of /i 7^ 1 should 
require extensive numerical studies that is out of our current interest. The variable Xn 
is decoupled from Yn corresponding merely to the linear chain, while Y^ represents the 
stretching motion that is our main interest. So, from now let us ignore the X^ part in 
further calculation. 

3. Mechanical statistics calculation 

Now let us calculate the mechanical statistics behavior of current system using transfer 
integral method. Having certain hamiltonian of an equilibrium system, one can extract 
some physical observables through statistical mechanics approach. The approach is 
particularly suitable to investigate the dynamics of DNA with thermal fluctuations since 
the exact solution of its equation of motion never exists [14] . Within the PB model, the 
DNA denaturation and its melting temperature have previously been studied in term 
of temperature [10, 12]. It was argued that the denaturation can be induced by energy 
localization due to nonlinear effects, and also should be highly influenced by external 
interactions like Vn in the present case. 

Using the hamiltonian in Eq. (2), one can consider the partition function to 
further calculate some thermodynamics variables. In [10, 14], the calculation has been 
performed using transfer integral method for V = 0. This paper follows the same 
procedure. In the canonical ensemble the partition function density is related to the 
Hamiltonian by Z oc exp{/3H). For the present case it reads, 

N N 

^dYn exp { -/ 



n=l n=0 



P^ K., ,2 



5^ 2M + 2 " ^'^-'^^ 
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+ f (e-^/''^" - 1)' + V{Y^^ I 5{Y^ - Yo) , (3) 

where /3 = 1/kT with k is the Bohzmann constant. It can further be decomposed into 
its momentum and coordinate spaces, Z = ZpZy with, 



N 

Y\ dP„ exp 

n=l 

TV 



'^^9M ) ' 



(4) 



Zy= fll dYn 6{Ym - lo)exp [^(F, - Y^-if 



n=0 

The Gaussian integration in Eq. (4) yields, 

N 

'27rM' 



Zp^ 



On the other hand, Eq. (5) reads. 



N 



N 



where. 



Zy= jt[ dr„ 5{Y^ - Yo) n e-'^®^^"'^"-) , 

n=0 n=l 



(5) 



(6) 



(7) 



(8) 



This can be solved using transfer integral method. 

Shifting the operator Y^-i — >■ Y^ and defining an eigenfunction 0(1^) that satisfies 

[14], 

J dy„_i e-^®(^-^"-)0,(y„_o = e-^^^0,(y„) , (9) 

Eq. (7) becomes. 



AT 

Zy= H dY^ 6{Yr, - Yo) e-^^^UYn) 

n=0 



(10) 



Since S{Yn — Yq) — X^i (^Jv)^j(^o) and taking the normalization J dY (f)*(f) — 1, one 
gets, 

Zy = J]e-^^^^. (11) 

i 

Let us find the eigenfunction (f){y) and the eigenvalue E of Eq.(9). First of all at 
continuum limit, Y^-i ~ x and F„ ~ y, Eq. (9) reads. 



J dx 0(x)exp I —/3 



2 (1/ - + f (e""/'" - 1)' + V{y) I = e-^^ 0(y) .(12) 
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Now assuming that the harmonic term is dominant, x — y -\- z can be expanded as a 
Taylor series in the power of z. The expansion up to the second order yields, 



• d0 IdV 

.^^^)+d^"+2d7 



ydze "2 

= exp|-/3 E-^{e-'^/^y-lf-rV{y) | 



(13) 



Again, this is just the Gaussian integration and the result is [15], 



0(y) + 



1 d20 

2^dy2 



= exp 



-/3 



-^(e-t^-l)VnZ/) 



0(Z/) 



0(2/) 



+ /9 



^(e-/^^-l)%F(l/) 



0(2/) 



This is nothing else than the Schrodinger like equation, 

1 a"0 U r „,/o,. 



(14) 



(15) 



where mo = and E = E +1/ (2(3) In [27r/(/3«;)]. The problem is therefore turned into 
finding the eigenvalue of Eq. (15). 

Eq. (15) then yields an eigen equation Hcf)^ — Em(f)m with Hamiltonian, 



H 



1 



D 



[e-'^/'y-l]' + V{y), 



(16) 



2mo dy2 2 

and m is an integer. For a constant force Vq, l^(y) = Voy, the transformation 
</) = exp{l3Voy)ip reproduces the Schrodinger equation with Morse potential [16]. 
Unfortunately, the equation can be solved only for few special cases. Generally, the 
solution can be obtained pcrturbatively by assuming that the potential V{y) is small 
enough compared to the Morse potential. Under this assumption, the eigen equation can 
be solved using time independent perturbation theory in standard quantum mechanics. 
However, it should be emphasized that the problem remains a classical one, and the 
quantum mechanics is borrowed only for technical reason since the calculation yields 
the Schrodinger equation. 

In quantum mechanics the perturbation is applied by expanding the Hamiltonian 
(H) in a form of H = Hq + eV with Hq is the hamiltonian in Eq. (16) without the 
potential, while e is a small parameter. According to the time-independent perturbation 
theory, the eigenvalue and its eigenfunction can be written as [17], 

+ ..., (17) 

(18) 



En 



+ eE 



(1) 



(2) 



My) = <t>^:Hy) + ^<l>il.Ky) + ^'€Ky) + -- 

The eigenvalue for Hq is determined by i?o0*^°^ = -E'o°V^°\ that is, 
1 d2(/.W D 



2mo dy^ 



+ _(e-"/2i/_ 1)^0(0) ^^(0)^(0)^ 



(19) 
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where E^^^ = £'(°) + l/(2^) ln[27r/(/3K)]. The equation is just the well known Schrodinger 
equation with Morse potential and the solution has been derived in [14]. Substituting 
y' = {l/2)ay, A = 2y/moD/a, S^^^ = Smo/a'^E^^^ and z ~ 2Aexp(— y'), one arrives at 
the eigenvalue problem for KD/a > 1/2 as [18, 17], 

(e-V - 2e-') <t>t' = . (20) 



dV ,2 



dy'2 

The eigenvalue and its eigenf unction are, 

2 



4°)=-(A-m-i) , (21) 



2^ 

(t>t^ =N^z'-e-^"Ll-{z), (22) 

where 6^ = A — m — 1/2 and = "m)- [r (m + 1) T (2A — m)]"^''^ with V is the Gamma 
function. L'^"^{z) = {z^'^^"'e^ /m\){d"^{e^^z™-^'^^"')/dz^) is the associated Laguerre 
polynomial or generalized Laguerre polynomials. 

Based on these transformation, using the standard procedure in quantum mechanics 
one obtains the solution up to the first order [17], 

K{y') 4>t\y') + ^4>t\y') . (24) 



where. 



41' = <W<<Ti.v')V(v'WS(v'), (25) 



(26) 

Finally the partition function of y is given by, 

= ^e-^^^'", (27) 

m 

with the eigenvalue, 

4. Thermal denaturation 

One of relevant parameters in the study of DNA denaturation is the mean stretching, 
, of the hydrogen bond. It is defined as, 

1 f ^ 

iVn) = ^ / dyndpn Hyne'"'' ■ (29) 

n=l 

The average value of hydrogen bond stretching can be calculated through. 
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At the limit of large A^, {y) is dominated by the ground state [14], that is 
{y) ~ I dy(0o(y))*y0o(y)- Up to the leading term of perturbation, i.e. (po{y) — 
0i°^(l/) + etf^(l/),Eq. (30) reads, 

/dy0r(?/)#r(?/) 

+ e J dycj>^^>{y)ycj>^^\y) + J dy 
The higher order of eigenvalue and eigenfunction are given by, 

- / d#r*(z/)l^(z/)0r(l/) , 



(31) 



(32) 



and. 



/dy0r\/(y)0r^(o) 



(0) 



0r(z/)+ 



/d^0|n%)0r 

c(0) _ -r{0) 



Then Eqs. (22)-(26) yield, 
(y) = ^(00) + e 



'(10) 



ciO) _ p(0) 
Cq 



'(20) 



c-(O) _ _f.(0) 
^0 ^1 



'(20) 



c-(O) _ c-(O) 
'-'0 '-'2 





^1 

(t>T* {y)y(t>t\y) 
dy (l>o^*{y)y(f>f\y) 

j dy4>?^\y)y<t>f{y) + --- +0{^) 



<t>T{y) + --- -(33) 



where, 

-^00 
-^20 

while the eigenvalues and eigenfunctions are given by, 

2 



j dyct>'S^\y)V{y)ct>f\y), 
j dy<t>f>{y)V{y)<j>f{y), 

dy^t^\y)V{y)ct>f{y), 



'(0) 



- A- 



1 



'(0) 



(34) 

(35) 
(36) 
(37) 

(38) 
(39) 
(40) 
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Here, 



(2A) 



A-1/2 



g-Aexp(-a/2 2/) g-a/2(A-l/2)j/ 



2(2A)^-i 

v/r(2)r(2A - 1) 

X [-Ae-"/2^ + (A-l)] , 
4(2A)^-5/2 

Vr(3)r(2A - 2) 

X 



(41) 



g-Aexp(-a/2 2/) g-a/2(A-3/2)j/ 



^_Aexp(-a/2;/) g-«/2(A-5/2)y 



A^ e""^ - A(2 A - 3) e'^/^ 2/ + 1 (2A - 3) (2A - 4) 



(42) 



(43) 



where A = 2\fiW / {aksT) . 

Considering the case of hnear potential V{y) — V^y, the wave function is given by, 



My) 



(0) 



(0) 



.(0) 



'(0) c(0) 



The variation of wave function in respect to temperature is depicted in Fig. 1 for the 
same parameters as [12], that is D = 0.04 eV, a = 4.45 A-^ e = 0.01, n = 0.06 eV/A^ 
and Vq = 15 pN. 

The ground state of eigenfunction for hnear potential gives the weighting factor for 
the calculation of DNA displacement (y) [10]. The figure shows that at low temperature, 
the peak is close to the minimum of Morse potential and decays exponentially as the 
displacement increases. Without the external potential, the eigenfunction is increasing 
drastically at the melting temperature, that is about 500K as already pointed out in [12]. 
On the other hand, the same behavior occurs at lower temperature when the external 
potential is applied to the system, i.e. at about 200K. 

Prom Eq. (34), the influence of various external potentials into denaturation 
processes is depicted in Fig. 2 using the same parameters as [12]. The figure describes 
the effects of hnear potential Vijj) = V^y, sinusoidal potential V{y) — yosin(ay) and 
Gaussian potential V{y) — T^e""^ with a is a constant. It should be noted that the 
result coincides to the previous work for zero external potential [12]. The figure shows 
that the critical temperature is approximately around 400K without external potential, 
and around 200K with constant external potential. Nevertheless the critical temperature 
becomes around 350K for periodic sinusoidal force with a = 2.5. The same behavior of 
(y) for non forcing condition has also been obtained in [19] using fiexible chain model 
with typical parameter D — 0.25ey, a — 2.8 A~^ and numerical calculation of the 
stochastic description of fiexible polymer models. 



5. Summary 



The effect of external potential to the Peyrard-Bishop DNA denaturation has been 
investigated for various types of potential : linear potential V{y) = Voy, sinusoidal 
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Figure 1. The eigenfunction without (left) and with (right) external linear potential 

as a function of the displacement in various temperatures : lOOK (solid line), 300K 
(dashed-line) and 500K (dotted-line) in the left; and lOOK (solid line), 150K (dashed- 
line) and 200K (dotted-line) in the right figures. 

potential V{y) = Vq sin(ay) and Gaussian potential Viy) = Voe~"^^ with a is a constant. 
The DNA denaturation and its melting temperature have further been studied in the 
framework of statistical mechanics approach. The calculation of partition function has 
been performed perturbatively using transfer integral method to obtain analytically the 
mean stretching of hydrogen bond. 

The present paper support the results of previous works. In a cell, DNA strands can 
be separated by applying certain external potential [20], or in chemical terms by enzymes 
whose interactions with DNA could make the strand separation thermodynamically 
favorable at ambient temperature [21]. It has been shown that two strands of double- 
stranded DNA can be separated by applying Vq ~ 15 pN force at room temperature. The 
model also predicts that the DNA overstretching force should be a decreasing function 
of temperature, or in another word the melting temperature should be a decreasing 
function of the applied force [20, 22]. 
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